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Abstract. We present 10—15 fim spectra of a sample of H ll regions, YSOs and evolved stars that show strong 
unidentified infrared emission features, obtained with the ISO/SWS spectrograph on-board ISO. These spectra 
reveal a plethora of emission features with bands at 11.0, 11.2, 12.0, 12.7, 13.5 and 14.2 fim. These features are 
observed to vary considerably in relative strength to each-other from source to source. In particular, the 10—15 
/xm spectra of the evolved stars are dominated by the 11.2 fim band while for H ll regions the 12.7 is typically as 
strong as the 11.2 band. Analysing the ISO data we find a good correlation between the 11.2 /xm band and 
the 3.3 /im band, and between the 12.7 fim and the 6.2 /im band. There is also a correlation between the ratio of 
the UIR bands to the total dust emission and the 12.7 over 11.2 /xm ratio. Bands in the 10—15 ^m spectra region 
are due to CH out— of— plane (OOP) bending modes of polycyclic aromatic hydrocarbons (PAHs). We summarise 
existing laboratory data and theoretical quantum chemical calculations of these modes for neutral and cationic 
PAHs. Due to mode coupling, the exact peak position of these bands depends on the number of adjacent CH 
groups and hence the observed interstellar 10—15 fim spectra can be used to determine the molecular structure 
of the interstellar PAHs emitting in the different regions. We conclude that evolved stars predominantly inject 
compact 100—200 C-atom PAHs into the ISM where they are subsequently processed, resulting in more open 
and uneven PAH structures. 

Key words. Circumstellar matter — Stars: pre-main sequence — H ll regions — ISM: molecules; — Planetary 
nebulae: general — Infrared: ISM: lines and bands 



1. Introduction 

Many sources show a rich set of emission features at 3.3, 
6.2, 7.7, 8.6 and 11.2 /im commonly called the unidenti- 
fied infrared (UIR) emission features. These bands domi- 
nate the IR spectra of a wide variety of objects, including 
H II regions, PNe, post-AGB objects, YSOs, the diffuse 
ISM, galaxies and galactic nuclei ( Cox fc Kessler 1999|) 



Laboratory spectroscopy of PAHs shows that, besides 
the well known UIR bands, PAHs exhibit many weaker 
emission bands. In particular, the region from 10 to 15 
/im has a rich spectrum due to the out— of— plane bend- 
ing vibrations (OOP) of aromatic H-atoms. The peak 
wavelength of these modes depends on the structure of 
the molecule; in particular on the number of neighbour- 
ing H-atoms per ring (e.g. Bellamy 1958 ; Hudgins & 



tion of their carriers has become an important problem 
in astrophysics. These UIR bands are generally attributed 
to IR fluorescence of small (~ 50 C-atom) polycyclic aro- 
matic hydrocarbons (PAHs) molecules. 



Under^jtanding these features and the origin and evolu- AUamandola 1999|). Here we present data in this region 



obtained with the Short Wavelength Spectrometer (SWS) 
(de Graauw et al. 1996) on-board the Infrared Space 



Send offprint requests to: S. Hony (hony@astro.uva.nl) 

* based on observations obtained with ISO, an ESA project 
with instruments funded by ESA Member states (especially the 
PI countries: France, Germany, the Netherlands and the United 
Kingdom) and with the participation of ISAS and NASA. 



Observatory (ISO) (Kessler et al. 1996). The sensitivity 
and medium resolving power of the instrument allows us 
to detect and resolve several weak features predicted by 
the PAH hypothesis and to determine the molecular struc- 
ture of the emitting PAHs and their evolution in space. 

In Sect. H, we present the observations of our sample of 
H II regions, YSOs, reflection nebulae (RNe) and evolved 
objects. The 10—15 /im regions of these sources are anal- 
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Fig. 1. Spectra of 3 sources that show features in the re- 
gion of interest. The dashed hues are the continua men- 
tioned in the text. 

ysed in Sect. |[ The spectral characteristics of PAHs in 
this wavelength range as measured in the laboratory and 
calculat ed b y quantu m c hemical theories are summarised 
in Sect. 4.1. In Sect. 4.2 the laboratory spectra are com- 
pared to the observed spectra. The molecular structures 
implied by the observed spectra are discussed in Sect. || 
while in Sect. ^ the origin and evolution of these molec- 
ular structures are examined. Finally in Sect. ^ our main 
results are summarised. 



2. Observations 

We include 29 high S/N spectra of 16 sources taken from 
the samples of H ii regions, YSOs, and post-AGB objects. 
For the trend analysis presented in Sect. 3.4 we have added 



another 15 sources in which the S/N is such that we cannot 
reliably measure the weakest features but we determine 
the intensities of the stronger UIR bands with reasonable 
accuracy. All spectra presented here have been obtained 
with the SWS instrument in the AOT01/AOT06 scan- 
ning mode at various speeds, with resolving power (A /AA) 
ranging from 500—1500 (see de Graauw et al. 1996| ), see 
Table |l] for details of the observations. 

The data were pro cessed using SWS inter active anal- 
ysis product; lA'^ (see ie Graauw et al. 1996 ) using cali- 
bration files and procedures equivalent to pipeline version 
7.0. If a source has been observed multiple times and these 
observations are of similar quality and of comparable flux- 
level these data are co-added after the pipeline reduction. 



Since the features we discuss here are fully resolved in 
all observing modes, we combine the data obtained in all 
different modes to maximise the S/N. Further data pro- 
cessing consisted of bad data removal, rebinning on a fixed 
resolution wavelength grid, removing fringes and splicing 
of the sub-bands to form a continuous spectrum. 

For all spectra the amount of shifting between sub- 
bands required falls well within the calibration uncertain- 
ties in the region of interest: 7 to 16 [ivn. Any jumps be- 
tween bands are due to either flux calibration and dark 
current uncertainties. The effect of dark current is most 
important in low flux cases while the flux calibration un- 
certainties will dominate in bright sources. Below 20 Jy 
we apply offsets to correct for dark current uncertainties. 
In these low signal cases the typical noise level in the dark 
current measurements of 1—2 Jy introduces offset uncer- 
tainties >5— 10 per cent dominating the flux calibration 
uncertainties. Above 20 Jy we apply scaling factors to cor- 
rect for flux calibration uncertainties. The splicing intro- 
duces little uncertainty in the measured strengths since 
most features fall completely within one ISO/SWS sub- 
band. An exception to this is the band strength of the 
12.7 /im feature. This feature is sensitive to the way band 
2C (7 to 12.5 /im) and 3A (12 to 16.5 /im) are combined. 
This introduces an extra uncertainty of the band strength 
of typically 20—30 per cent for the weakest features. 

Some SWS data, especially in band 3A, are affected 
by fringes. We have corrected for fringes in those sources 
where they occur, using the aarfringe tool of lA'^ on the 
rebinned spectrum. In the method we apply fringes are 
fitted with sine functions with periods in the range where 
fringes are known to occur and divided out. Note that the 
features we study here are much broader than any of fringe 
periods, therefore the intensities we measure are not di- 
rectly affected by the fringes. However in some cases after 
fringe removal the continuum is more easily determined. 

SWS spectra of many sources, including stars en- 
shrouded in both carbon-rich and oxygen-rich dust and 
sources without any circumstellar material show very weak 
structure around 13.5 and 14.2 /im at the 3 to 4 per cent 
level relative to the continuum possibly due to residual in- 
strumental response. The emission features discussed here 
are all stronger than this with a maximum of 85 per cent of 
the continuum in the reflection nebula NGC 7023. Near 
11.03 /im there is a residual instrumental feature which 
coincides with the weak 11.0 /im feature that we observe 
in our spectra. We have included the effect of this feature 
in the uncertainty on the intensities in Table |2[ 

Many sources in this sample have strong narrow emis- 
sion lines in their spectrum, in particular the strong [Ne ll] 
line at 12.81 /im is perched on top of the 12.7 /im UIR 
band. This line and the UIR band are easily separated 
at the resolution of the SWS instrument. We remove the 
contribution from this line by fitting a Gaussian profile 
to the line and subtracting that profile prior to rebinning. 
The spectrum of NGC 7027 has a very strong [Ne V] emis- 
sion line at 14.32 /im. We have removed the part of the 
spectrum which contains this line. 
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Table 1. Source list. Observational details of the sources used in this study. 
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Object 


Obs." 


Q 


5 


TDT" 


Sp.Type 


Go 


Obj.Type 




M ode 


(J2000) 


(J2000) 




(K) 


(1.610 




AFGL 437 


01(2) 


03 07 23.68 


+58 30 50.62 


86300810 


08.5 


1E5 


Star forming region 


IRAS 03260+3111 


01(3) 


03 29 10.37 


+312158.28 


65902719 


B9 


2E4 


Herbig AeBe 


NGC 2023 


01(3) 


05 4138.30 


-0216 32.59 


65602309 


B1.5V 


3E2 


Refi. Nebula 


HD 44179 


01(4) 


06 19 58.20 


-10 3815.22 


70201801 


B8V 


5E6 


post-AGB 


IRAS 07027-7934 


01(2) 


06 59 26.30 


-79 38 48.01 


73501035 


WCIO 


2E7 


PN 


HD 97048 


01(4) 


1108 04.61 


-77 3916.88 


61801318 


AO 


2E4 


Herbig AeBe 


IRAS 12405-6238 


01(3) 


12 43 31.93 


-62 55 11.39 


29400410 


09.5 


3E5 


H II 


HEN 2-II3I" 

- 


01(1) 


14 59 53.49 


-5418 07.70 


07903307 


WCIO 


6E4 


PN 




01(2) 


14 59 53.49 


-5418 07.70 


43400768 








IRAS 15384-5348 


01(2) 


15 4217.16 


-53 58 31.51 


29900661 




1E4 


H II 


CD-42 11721(off) 


01(2) 


16 59 05.82 


-42 42 14.80 


28900461 


BO 


— 


Herbig AeBe (off pointing) 


CD-42 11721"^ 










BO 


- 


Herbig AeBe 


- 


01(2) 


16 59 06.82 


-42 42 07.60 


08402527 








- 


01(2) 


16 59 06.80 


-42 42 07.99 


64701904 








IRAS 17047-5650^ 










WCIO 


5E6 


PN 


- 


01(3) 


17 09 00.91 


-56 54 47.20 


13602083 








- 


01(1) 


17 09 00.91 


-56 54 48.10 


27301339 








HB 5 


01(3) 


174756.11 


-29 59 39.70 


49400104 


120 OOOt 


- 


PN 


NGC 6537 


01(3) 


18 0513.14 


-19 50 34.51 


70300475 


AO 


- 


PN 


GGD 27-ILLt 










Bl 


3E6 


Star forming region 


- 


01(2) 


18 1912.04 


-20 47 30.98 


14802136 








- 


01(2) 


18 1912.03 


-20 47 30.59 


14900323 








IRAS 18240-0244 


01(1) 


18 26 40.00 


-02 42 56.99 


14900804 


WC8 


- 


PN 


IRAS 18317-0757 


01(2) 


18 34 24.94 


-07 54 47.92 


47801040 


08 


- 


H II 


IRAS 18416-0420 


01(2) 


18 4415.19 


-0417 56.40 


13402168 


05.5 


- 


H II 


IRAS 18502-K0051 


01(2) 


18 52 50.21 


+00 55 27.59 


15201645 


07 




H II 


TY CRA^ 










B9 


6E3 


Herbig AeBe 


- 


01(3) 


19 0140.71 


-36 52 32.48 


33400603 








— 


01(1) 


19 0140.70 


-36 52 32.59 


34801419 








- 


01(3) 


19 0140.71 


-36 52 32.48 


71502003 








BD +30 3639 


01(3) 


19 3445.20 


+30 30 58.79 


86500540 


WC9 


1E5 


PN 


IRAS 19442+2427 


01(2) 


19 46 20.09 


+24 35 29.40 


15000444 


07 


6E6 


H II 


BD+40 4124 


01(3) 


20 20 28.31 


+412151.41 


35500693 


B2 


1E4 


Herbig AeBe 


S106 IRS4 


01(2) 


20 27 26.68 


+37 22 47.89 


33504295 


08 


1E5 


H II 


NGC 7023^ 










B3 


5E2 


Refi. Nebula 


- 

- 


01(4) 
01(2) 


210131.90 
210130.40 


+68 10 22.12 
+6810 22.12 


20700801 
48101804 








NGC 7027+ 










200 OOOo 


2E5 


PN 


- 


01(4) 


210701.71 


+42 14 09.10 


02401183 








- 


01(1) 


210701.70 


+42 14 09.10 


23001356 








- 


01(2) 


210701.70 


+42 14 09.10 


23001357 








— 


01(3) 


210701.70 


+42 14 09.10 


23001358 










06 


210701.50 


+4214 10.00 


33800505 










01(4) 


210701.63 


+421410.28 


55800537 








IRAS 21190+5140 


01(2) 


2120 44.85 


+5153 26.59 


15901853 




2E5 


H II 


IRAS 21282+5050''^ 










09 


1E5 


PN 




01(2) 


2129 58.42 


+5103 59.80 


05602477 










01(3) 


2129 58.42 


+5103 59.80 


15901777 










01(2) 


2129 58.42 


+5103 59.80 


36801940 








IRAS 22308+5812 


01(2) 


22 32 45.95 


+58 28 21.00 


17701258 


07.5 


3E3 


H II 


IRAS 23030+5958 


01(2) 


23 05 10.57 


+60 14 40.60 


22000961 


06.5 


7E3 


H II 


IRAS 23133+6050 


01(2) 


23 15 31.44 


+6107 08.51 


22001506 


09.5 


7E5 


H II 



° SWS observing mode used (see ^^^^S^^^^L 1996| ). Numbers in brackets correspond to the scanning speed. 
'TDT number which uniquely identifies each ISO observation. 

+ These spectra have been ob tained by co-adding the separate SWS spectra also list ed in the table, see text. 
^Effective temperature from Gesicki & Zijlstra (2000). ^Effective temperature from Latter et al. (2000). 
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We also include in Table |i| the spectral type of the il- 
luminating source and an estimate of the flux density at 
the location where the PAH emission originates from in 
units of the average interstellar UV field ( Habing 1968| ). 
We have derived these estimates from the observed IR 
flux (/ir) and the angular size of the PAH emission re- 
gion (Wolfire et al. 1989). This estimate is based on the 
assumption that all the UV light is absorbed in a spherical 
shell with the angular size of the object and re-emitted in 
the IR. The flux density at the shell is given by: 



Go 



1.610- 



1 pc 
1 AU 



02 



1.110^^61"^ 



(1) 



where Go is the UV flux density in 1.6 10~^ Watt/m^ and 
is the angular diameter of the object in arcseconds. We 
have used for the size of the H ii regions the measured ra- 
dio sizes. This is reasonable since the PAHs are expected 
to be destroyed inside the H ii region. For the other ob- 
jects we estimate the size of the PAH emission region from 
ISOCAM data ( |Cesarsky et al. 1996 ) when available. 



3. Results 

3.1. Overview 

In Fig. |l| we present spectra of three typical sources to 
illustrate the spectral detail present and the variety of un- 
derlying continua. Even more so than previous IRAS /LRS 
and ground-based studies suggested ( Cohen et al. 1985| ; 
Roche ct al. 1989; Witteborn et al. 1989), the complete 



11 — 15 iim spectra reveal an extremely rich collection of 
emission features with bands at 10.6, 11.0, 11.23, 12.0, 
12.7, 13.5, and 14.2 /im. These features are perched on 
top of an emission plateau of variable strength, which ex- 
tends across the entire region. 

3.2. Continuum 

In order to compare the intrinsic strength of the UIR 
bands in the different sources, we subtract a continuum, 
splined through points from 9—10.5 and 14.5—15.5 /im and 
through points near 11.8 and 13.1 /im. Since GGD 27 and 
IRAS 18416 have broad solid CO2 absorption features be- 
yond 15 fj,m, we take 14.7 /im to be the continuum rather 
than extending to 15.5 fim. Typical examples of the con- 
tinua are included in Fig. |l|. Besides the chosen points, the 
continuum also runs through the local minimum at 10.9, 
12.2 and 14.0 /tm with the exception of the Red Rectangle. 
This object has excess emission around these wavelengths, 
which could be related to the presence of crystalline sili- 



cates in the vicinity of HD 44179 ( [Waters et al. 19981 ). For 
any reasonable choice of continuum points, the underlying 
plateau ranges from about 10 to 14 /im, peaking at about 
12 /im. At first sight, the plateau of GGD 27 seems to dif- 
fer with an onset at about 11 /xm and peaking around 13 
/im. However, this likely reflects the effects of foreground 
sihcate absorption. The PN IRAS 21282-^5050 does show 
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Fig. 2. An overview of the observed features. All spectra 
have been continuum subtracted and are scaled to have 
the same integrated intensity in the 12.7 /im feature. The 
sources are ordered according to their 11.2/12.7 /im band 
strength ratio (bottom to top). The ratio of the 11.2 /tm 
to the 12.7 /im feature spans a full decade. 
^ Sources with broad solid CO2 absorption beyond 15 /im. 



a plateau that differs, peaking at 12 /tm but extending all 
the way to 17 /tm. 

3.3. Emission bands 

We show the continuum subtracted spectra of the 16 ob- 
jects with the highest S/N in Fig. ^ Perusal of the spec- 
tra reveals a plethora of weaker features. NGC 7027, HD 
44179 and IRAS 21282 have a very weak, broad feature 
near 10.6 /im. Almost all these sources show a feature at 
11.0 /im, except for GGD 27, where the feature is possibly 
present but only at the la level. The weak 12.0 /tm band 
is detected in 10 out of 16 sources. The 13—15 /tm range 
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Table 2. Peak position and strength of the features observed in the 10—15 fim for the sources shown in Fig. 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10)^ 

/l2.0 * 


(11) 


Source 


AAii.at 
[lO^^^im] 






Ac, 13. 5 


/l3.5 


Ac, 14. 2 

[fj.m] 


/l4.2 


Jll.O 


r bo 
-t plateau 


IRAS 03260 


0.9(0.1) 


15.8(0.4) 


6.4(0.4) 


13.57(2) 


1.7(0.1) 


14.19(2) 


0.9(0.1) 


0.8(0.1) 


n 


22 


HD 44179 


12.6(0.5) 


118.1(6.4) 


26.3(1.4) 


13.61(2) 


13.4(0.5) 


14.21(2) 


8.6(1.3) 


4.8(2.3) 


5 


112 


HEN 2-113 


2.7(0.3) 


17.0(1.0) 


5.6(1.0) 


— 


— 


14.26(5) 


4.7(1.8) 


0.5(0.5) 


d 


30 


IRAS 15384 


0.0(0.6) 


14.5(0.8) 


13.5(0.6) 


13.52(2) 


1.5(0.3) 


14.30(2) 


2.7(0.1) 


0.6(0.1) 


2 


66 


CD -42(off) 


-6.5(0.2) 


22.2(0.6) 


11.3(0.5) 


13.55(2) 


1.7(0.5) 


14.22(2) 


0.8(0.2) 


1.0(0.1) 


2 


25 


CD -42 


-8.2(0.8) 


31.7(1.5) 


16.7(1.0) 


13.53(2) 


1.6(0.6) 


14.22(2) 


1.6(0.3) 


1.9(0.3) 


1 


58 


IRAS 17047 


11.8(0.4) 


25.1(1.0) 


11.6(1.0) 


13.50(3) 


1.8(1.3) 


14.23(3) 


2.2(1.0) 


1.3(0.8) 


n 


145 


GGD 27 


2.6(0.6) 


8.9(0.6) 


5.4(0.5) 


13.57(2) 


2.1(0.2) 


14.27(4) 


1.5(0.9) 


1.4(0.1) 


d 


3 


IRAS 18317 


6.5(0.1) 


10.5(0.2) 


15.6(0.3) 


13.44(2) 


2.0(0.3) 


14.22(2) 


2.2(0.2) 


0.2(0.1) 


1 


83 


IRAS 18416 


-2.3(0.3) 


12.8(0.9) 


9.8(0.5) 


13.55(2) 


2.4(0.4) 


14.18(3) 


3.1(1.2) 


1.2(0.3) 


n 


52 


TY CRA 


-5.0(0.3) 


15.5(0.1) 


4.7(0.2) 


13.54(2) 


1.3(0.3) 


14.21(2) 


0.5(0.1) 


0.8(0.1) 


n 


13 


S 106 


-4.1(0.1) 


19.7(0.9) 


15.6(1.0) 


13.56(2) 


2.7(1.2) 


14.19(2) 


3.6(0.2) 


1.1(0.3) 


2 


61 


NGC 7023 


-3.9(0.2) 


9.6(0.8) 


4.3(0.3) 


13.50(2) 


1.1(0.2) 


14.21(2) 


0.7(0.3) 


0.5(0.1) 


n 


11 


NGC 7027 


2.8(0.2) 


142.7(5.4) 


35.9(2.5) 


13.52(2) 


9.7(1.5) 


14.26(2) 


11.1(1.5) 


4.0(1.7) 


6 


451 


IRAS 21282 


4.2(1.4) 


20.4(0.7) 


6.2(0.2) 


13.40(2) 


1.5(0.4) 


14.22(2) 


1.7(0.6) 


0.6(0.3) 


y 


104 


IRAS 22308 


-3.2(0.5) 


8.9(0.4) 


5.8(0.6) 


13.55(4) 


1.3(0.6) 


14.23(2) 


2.1(0.1) 


0.7(0.1) 


n 


9 



The intensities /11.2, /12.7, Iii.o, I12.0 and /piatcau are the integrated fluxes of the features after continuum subtraction (columns 
3,4,9,10 and 11). The position of the 11.2 and 12.7 jj,m features are determined by fitting the template profile to the spectra in 
which we allow for a wavelength shift and a scaling of the template (column 2). The properties of the features around 13.5 and 
14.2 are determined by fitting Gaussian profiles to the data (columns 5—8). Numbers in parentheses are uncertainties. 
"Central wavelength of the fitted Gaussian profile, uncertainties are given in parentheses in units of W~^fim. 
''Intensities and uncertainty units are 10~^''W/m^ 

^Shift of the 11.2 ^m feature relative to the template 11.2 /xm profile. The template profile peaks at 11.229(0.001) fim. 

^No significant shifts of the 12.7 /im profile are observed with respect to the mean 12.7 /im profile. The mean profile ranges 

from 12.2—12.95 pm with the peak position at 12.804(0.005) /im. Typical uncertainty in the position determination 0.006 /im. 

* Typical error on 7i2.o is 0.5-10-'''W/m^ d means detected but not measured; n means not detected. 

^Typical error on /piatcau is 10T0~^*W/m^. 



contains the newly discovered weak 13.5 and 14.2 fim fea- 
tures (Fig. H). Structure near 13.5 is present in all sources 
except HEN 2-113, however for IRAS 21282 and IRAS 
18317 the band is replaced by a feature that peaks at a 
shorter wavelength. A feature near 14.2 fim is found in all 
sources although the actual peak position varies consider- 
ably ( cf., IRAS 03260 and IRAS 15384). This variation in 



peak p jsition of these two weak bands is quite in contrast 
with that of the 11.2 and 12.7 fim bands which almost in- 
variably peak at about the same wavelength (see below). 
A priori, it is not given that the spectral structure near 
14 jim is actually related in the different sources. Possibly, 
rather than one molecular vibration with a varying peak 
position, at this weak level of emission, we are probing 
different bands whose relative strengths reflect the condi- 
tions in the different sources. 



In Fig. ^ we show the continuum subtracted spectra 
after normalising to the integrated strength of the 12.7 
/xm feature. The sources are ordered according to the 
strength of the 11.2 /xm feature relative to the 12.7 fim 
band. Relative to the 12.7 /im band, the sources with the 
weakest 11.2 /xm feature are the H 11 regions (at the bot- 
tom of Fig. ^ , while the evolved stars show the strongest 
11.2 /im band. 



The spectral characteristics of the features are sum- 
marised in Table ^. Note that the uncertainties quoted in 
the Table reflect the noise level and the freedom in draw- 
ing the continuum within the methodology used to mea- 
sure these bands. Other ways of decomposing the broad, 



results (e.g. 


Boulanger et al. 1998 




Uchida et al. 


Verstraete et al. 2001 


). However these differences 



are systematic differences and do not affect the source-to- 
source variations we observe. The intensities of the 11.2 
and 12.7 fim features are obtained by direct integration 
above the chosen continuum. We measure the peak posi- 
tion of the 11.2 and 12.7 fim bands by fitting them with 
template spectra of these features. The template spectrum 
for the 11.2(12.7) fim feature is constructed by adding 
the continuum subtracted spectra with each 11.2(12.7) 
/im feature normalised to have the same integrated in- 
tensity. This way each source has equal contribution to 
the template spectrum. We use a -minimisation rou- 
tine to fit the template to the sources, allowing for both a 
wavelength shift and scaling in strength. The shifts that 
we determine for the 11.2 /xm band are very small except 
for HD 44179 and IRAS 17047 where this band is much 
broader than the template spectrum (cf.. Table |^, see also 
Peelers et al. (2001, in prep.)). Although there are differ- 
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Fig. 3. An overview of the observed features near 13.5 
and 14.2 fim. ^RB 44179 and NGC 7027 have been 
scaled by a factor of 0.3. 



ences between the detailed profiles of the 12.7 /im band 
we detect no significant shift of the band as a whole. For 
the weak features near 13.5 and 14.2 /Ltm, the parameters 
have been determined through fitting of Gaussian profiles. 
We adopted a local linear continuum for the very weak 
11.0 /im feature because of the severe blending of this 
band with the 11.2 fim band. The weak 12.0 /zm band is 
close to both the 11.2 and the 12.7 iim band. For only a 
few sources we measure the intensity of this band, for the 
other sources we refrained from detailed analysis. However 
Table |2| does note whether we detect this band. 

The profile of the 11.2 /im feature is asymmetric with 
a sharp blue rise and a more gradual decline to longe r 
wavelengths (|Roche et al. 1989| ; [Witteborn et al. 1989|) . 
This will be discussed in more detail for this sample by 
Peeters et al. (2001, in prep.). The 12.7 /im band is also 



asymmetric but in the opposite way with a slow blue rise 
and a sharp red decline between 12.8 and 12.9 /im. Because 
of their intrinsic weakness, the profiles of the 10.6, 11.0, 
12.0, 13.5, and 14.2 /tm features in the individual sources 
are not well determined however in the averaged spectrum, 
these features appear symmetric (cf., Fig. |^. 



3.4. Trends 

Given the large range of band strength ratios in the region 
of interest and the smooth variations in this ratio, it is of 
interest to investigate the correlation of the strength of 
all UIR emission bands with these variations. As a class 
all the planetary nebulae with a Wolf-Rayet central star 
(HEN 2-113, BD+30, IRAS17047, IRAS 07027 and IRAS 
18240) as well as the post-AGB star HD 44179 show a 
distinctly difi^erent UIR spectrum characterised by a shift 
in of the 6.2 and 7.7 /xm bands towards 6.3 and 8 /im 
respectively (Peeters et al., 2001, in prep.). Because of 
these spectral differences we have not included them in the 
following trend analysis in which they also separate out 
as a peculiar subgroup. These objects will be discussed in 
detail in a forthcoming paper (Hony et al., 2001, in prep.). 

Because here we want to study variations in the relative 
strength of the UIR bands to each other, not differences in 
absolute intensities differences due to intrinsic luminosity 
and distance of the source, we use 3-feature intensity ratio 
correlations. Although we observe variations in all ratios, 
we find only three that correlate and these are shown in 
Figures 16[ First, we find that the CH stretch mode at 
3.3 /(m correlates with the 11.2 /tm band (cf. Fig. ^). Note 
that the slope of the trend is roughly 1, which means that 
the I11.2/I3.3 is on the average constant at a value of 3—4. 

Second, the 12.7 /tm band correlates with the CC 
stretch mode at 6.2 /im, albeit with more scatter (cf. 
Fig. H) than the 11.2/3.3 ratio. The 7—9 /tm complex is 
well correlated with the 6.2 /jm band and similar plots 
can be made with these interchanged. Inspection of Fig. || 
and Table ^ shows that there is some indication for both 
the /13.5//11.2 and the IU.2/I11.2 /tm band to be higher in 
H II regions, however only about half the sources have such 
high S/N that these intensities can be reliably measured 
and this trend is not statistically significant. 

Lastly, we show in Fig. ^ the correlation between the 
ratio of the fiux emitted in the PAH bands relative to 
the total flux emitted in the IR (/ir) and the chang- 
ing /12.7/A1.2 ratio. We measure the /ir by integrating 
the SWS data and Long Wavelength Spectrometer (LWS) 
data if available. For those sources without LWS data we 
use a blackbody fitted to IRAS measurements in the wave- 
length region from 45—200 /tm. We do not apply correc- 
tions for aperture differences between the instruments. We 
estimate an uncertainty of 15 per cent on the /ir. Again 
different classes of objects occupy different parts in this 
diagram. 
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Fig. 4. Bands strength ratios as derived from the SWS 
spectra. Hexagons are H ii regions, stars intermediate 
mass star forming regions, squares RNe and triangles are 
PNe. 
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Fig. 6. The ratio between the flux emitted in the PAH 
bands over the total amount of IR radiation against the 
-fi2. 7/^^11.2 ratio. Plotting symbols are the same as in 
Fig. |. 



3.5. Comparison with other studies 

It is interesting to compare the result we obtained from 
the SWS spectra with results by other authors. Stu dying 
many locati ons in the diffuse interstellar medium Chan 



et al. (2000|) find that the relative bands strengths of the 
6.2, 7.7 and the 11.2 /im features do not vary system- 
atically over a wide range of intensities of the incident 
radiation field. Their observed ratios of 111.2/^6.2 cluster 
around 0.8, the value we observe in the RNe and the YSOs. 



Study ing a few RNe at various locations Uchida et al 
(2OOOD with ISO/CAM (ICesarsky et al. 1996| ) find only 



-0.6 -0.4 -0.2 

10glo(ll2.7/lll.2) 
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Fig. 5. Bands strength ratios as derived from the SWS 
spectra. Plotting symbols are the same as in Fig. 0. 



We also checked for correlations between band strength 
ratios and the flux density; Gq, however we do not detect 
any correlations. 

We emphasise that, while all UIR bands show a loose 
corr elation in the absolute intensity (see also Cohen et al 



small variations in band strength ratios. The /n. 2/^6. 2 ra- 
tios they derive with the method which is most similar 
to ours (method 2) are like those we find for the RNe. 
Their /12.7//11.2 ratios are however systematically lower 
that ours. This is due to the lower spectral resolution of 
the CAM spectra which results in blending and smearing 
of the weakish band at 12.7 fim. The SWS spectra are 
not affected by smearing since the features are fully re- 
solved. The results of these authors, like the results we 
present here, demonstrate that similar types of sources 
show similar PAH band strength ratios. Our results also 
demonstrate how different classes of objects show system- 
atic differences in their PAH spectra. 

Observations of the starforming region M17 have sug- 
gested that th e 13.5 fira band is corre lated with the mid- 
IR continuum Verstraete et al. (1996 ). Such a correlation 



1986^ ,989 ), these three are the only tight correlations 



present in this sample. 



is important to establish since it might yield information 
on the size of the carriers of the 13.5 /xm band. We have 
examined whether a correlation is also present in the sam- 
ple of sources we study here. We have therefore compared 
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Table 3. Wavelength region limits and the integrated ab- 
sorption cross-sections for the CH out— of— plane bending 
modes. 





Alow^ 






4 * 

^neutral 


A ■ * 

^cation 






[/xm] 


[km/mol] 


[km/mol] 


[km/mol] 


Solo 


10.6 


11.4 


24.8(13.5) 


25.7(14.2) 


24.1(12.9) 


Duo 


11.35 


12.8 


4(2.5) 


4.4(2.4) 


3.7(2.5) 


Trio 


12.5 


13.3 


9.6(5.9) 


10.1(5.3) 


9.0(6.5) 


Quartet 


13.0 


13.9 


12.0(4.8) 


11.5(5.5) 


12.6(3.9) 



Summary of the laboratory results on CH out— of— plane bend- 
ing modes for solo, duo, trio and quartet hydrogens on matrix 
isolated neutral polycy clic aromatic hydroca rbons and their 
cations. (Adapted from Hudgins et al. (2000b)). 
^Aiow is the lower limit of the region in fim. \up is the upper 
limit of the region in fim. 

'^The cross-section values for the solo, trio, and quartet modes 
per hydrogen are the averages over the spectra in the database. 
However, the A values for the duo mode per hydrogen decreases 
rapidly with size and settles to slightly less than 4 km/mole 
for PAHs with more than 24 carbon atoms. This value is more 
appropriate to use in determining the edge structures of PAHs 
that dominate emission in this wavelength region. 
*Total average cross-sections over both neutrals and cations in 
the database. 



the 15—16 /im continuum with the strength of the 13.5 /im 
band. We find strong variations, by a factor of ~100 in the 
strength of the mid-IR continuum relative to the 13.5 /im 
band. These strong variations are not surprising consider- 
ing the fact that we look at very diverse regions with large 
differences in dust composition and temperature distribu- 
tions. For example the evolved object HD44179 has con- 
tributions from crystalline silicates around 15 /im (Waters 
et al. l |998| ) . Many H ii regions have a strongly rising con- 
tinuum due to warm dust. However we would also like to 
point to the two observations of CD-42 11721 where the 
13.5 band is equally strong but the dusty continuum is 
missing in the off-pointed observation. These two observa- 
tions show the 15 micron continuum and the 13.5 micron 
feature to be decoupled even within the same object. 

4. The CH out— of— plane bending modes 

4.1. Laboratory spectroscopy of the OOP modes 

Chemists have long recognised the diagnostic value of the 
aromatic CH out— of— plane bending features in the 11 
to 15 fj,m spectral region for the classification of the aro- 
matic ring edge structures present in a particular sam- 
ple (e.g. Bellamy 1958). Specifically, the positions of the 
bands in this spectral region reflect the number of adjacent 
CH groups on the peripheral rings of the PAH structure 



(Bellamy 1958; AUamandola et al. 1985; Cohen et al. 1985; 



Leger et al. 1989|; [Roche et al. 1989|; |Wittcborn et al. 1989| ; 
AUamandola et al. 1999; Hudgins fc AUamandola 1999|) . 
Traditionally, aromatic rings carrying CH groups which 
have no neighbouring CH groups (termed "non- adjacent" 
or "solo" CH groups) show IR activity between 11.1 and 
11.6 /im. Likewise, activity between 11.6 and 12.5 /im is 
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Fig. 7. A comparison of the average interstellar spec- 
trum (top) with the ranges for the out— of— plane bending 
modes (bottom). The average spectrum was obtained by 
co-adding the continuum subtracted spectra after normal- 
isation to the 12.7 /im band strength. The boxes indicate 
the wavelength regions associated with the out— of— plane 
bending vibrations for different types of adjacent hydro- 
gen atoms determined from matrix isolated spectroscopy 
of neutral and cationic PAHs (see Hudgins et al. (2000b| ) 



for details). In this comparison it should be kept in mind 
that the emission process leads to a small (~ 0.1 /im) 
wavelength redshift in the peak position. 



indicative of two adjacent CH groups ("doubly-adjacent" 
or "duet" CH's) on the periphery of the PAH. Three adja- 
cent CH groups ("triply-adjacent" or "trio" CH's) are in- 
dicated by activity in the 12.4 to 13.3 /im region, and four 
adjacent CH groups ("quadruply- adjacent" or "quartet" 
CH's) by activity between 13 and 13.6 /im. Five adjacent 
CH groups ("quintuply-adjacent" or quintet CH's) are in- 
dicated by features falling in the 13 to 13.7 /im range. 
Trios and quintets also show a weak CCC bending mode 
in the 14—14.5 range. Other CCC bending modes occur 
in the 15—20 /im range and have been discussed in the 
astrophysical context by [Van Kerckhoven et al. (2000[ ) 



and Moutou et al. (2000). Over the years the reliability 
of this region to yield insight into the molecular structure 
and ring sidegroup placement on aromatic samples has 
been verified again and again (see Hudgins & AUamandola 
(1999D and references therein). However, most of these 



chemist's guidelines were based on studies of small PAHs 
where varying patterns of sidegroup substitution were 
employed to achieve different degrees of CH adjacency. 
Furthermore, these chemist's 'rules-of-thumb' are based 
on spectroscopic studies of aromatic molecules in solution 
or solid mixtures, environments quite different from that 
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of the emitting aromatic species giving rise to the interstel- 
lar features presented here. To obtain infrared data more 
relevant to the interstellar situation on larger molecules 
with the spectral details in this region determined by 
PAH structure rather than side-group substitution pat- 
tern, several groups have undertaken new spectroscopic 
studies of PAHs carried out under more appropriate con- 



not modified substantially compared to Bellamy (1958), 
ionisation causes some important changes in region bound- 
aries. These data expand on the initial report that the 
PAH cation solo hydrogen position is substantially blue 
shifted with respect to the wavelength for its neutral coun- 
terpart while the domains indicative of the other types 



of hydrogen are less affected by ionisation (Hudgins & 



ditions (e.g. see Szczepanski et al. (1995) and references Allamandola 1999). Considering these modified domains 



therein; Hudgins et al. (2000a) and references). Thanks to 
this effort, the mid-IR spectra of a few gas phase and many 
matrix isolated neutral and ionised PAHs are now avail- 



and taking into account the roughly 0.1 fim redshift in 
the peak position for PAHs emitting at temperatures of 



able. I i the Astrochemistry Laboratory at NASA Ames 1992 ; Colangeli et al. 1992| ; Joblin et al. 1995 ; Cook 



500 -1000 K (iFhckingcr et al. 1991|; [Brenner fc Barker] 



this w ork has been expanded considerably, including an 
extensive set of theoretical calculations, aimed to specifi- 
cally address the questions raised by these new ISO spec- 
tra. Since a detailed presentation and analysis of the 1 1 to 
15 iim region of the expanded dataset of matrix isolated 



Saykally 199^ ) allows us to draw the following conclusions. 



PAH spectra will be published separately ( Hudgins et al 
2000b^, only the salient points are summarised here and 
used in the analysis presented below. The IR spectra of 
matrix isolated PAHs compare favourably with the avail- 
able spectra of gas phase PAHs, for both neutral PAHs 
and cationic PAHs (see also Piest et al. 1999) and val- 
idates the use of the matrix isolation method to obtain 
astronomically relevant data. 

There are two points that emerge from an analysis 
of the laboratory database that are of particular impor- 
tance to the observational data presented here. The first 
involves the effect of ionisation on the characteristic wave- 
length regions of the various CH adjacency classes. The 
second is the intrinsic integrated absorption strengths (A 
values) which are derived for the various adjacency classes. 
Together these results provide the tools to not only qual- 
itatively infer the sorts of PAH edge structures present, 
but also quantitatively determine their relative amounts. 
As shown below, this allows one to place stringent con- 
straints on the emitting interstellar PAH family. 

4.2. OOP modes in the interstellar spectrum 

Fig. and Table | summarise the key points presented in 
Hudgins et al. (2000b ) that are applicable to this work. 
Fig. |7| schematically compares the average UIR spectrum 
with the wavelength regions associated with different hy- 
drogen type for neutral and ionised isolated PAHs, while 
Table ^ lists the specific wavelength limits and the in- 
tegrated band strengths per CH group as a function of 
hydrogen adjacency for all PAHs in the NASA Ames 
database. The regions indicated for the neutral PAHs dif- 
fer slightly from those indicated by Bellamy (1958| ). We 
deem our results more astrophysically representative be- 
cause of the larger set of molecules studied and because 
no substitution with strongly electro-negative groups were 
involved. Moreover our data were measured on isolated 
PAHs in inert matrices, rather than in solid mixtures. 

Perusal of Fig. and the wavelength limits listed in 
Table ^ shows that, while the ranges for neutral PAHs are 



— The broad, weak interstellar emission feature between 
10.6 and 10.7 fiia and the stronger distinct interstellar 
band peaking near 11.0 /im fall in the region unam- 
biguously attributable to PAH cations. 

— Adjusting for the 0.1 /im redshift, the bulk of the 11.2 
/im interstellar band falls squarely within the solo re- 
gion for the neutral PAHs and also, at the long wave- 
length end of the range for solo, cationic aromatic CH 
bonds. 

— The domains indicated in Fig. |^ show that regardless 
of the region definitions used, there can be little doubt 
that the weak interstellar 12 /tm band arises from duo 
modes. 

— Interestingly, the blueward skewed profile of the mod- 
erately strong band at 12.7 /im seems to fall better 
within the envelope for trio hydrogens in ionised PAHs 
and does not agree well with the envelope for trio 
modes in neutral PAHs. An origin in duo modes of 
neutral PAHs is also consistent (Fig. 0). 

— The weak 13.5 /im feature falls in the quartet domain. 
However, this also overlaps the lower of the two do- 
mains characteristic of the quintet region, and so, al- 
though we consider this highly unlikely, it is possible 
that quintet types of PAH structure could contribute 
to this band as well. 

The large laboratory database and the high quality of the 
ISO interstellar spectra allows us to compare the distri- 
bution of peak positions of the individual modes with 
the bands observed in the interstellar spectrum. From 
the measured peak positions and cross-sections from the 
Hudgins database we construct a synthetic spectrum of 
a mixture of PAHs for comparison with the observed 
interstellar UIR bands. We construct such a spectrum 
on a mode-per-mode basis. Per mode, we take for each 
molecule which shows this mode the measured peak posi- 
tion and cross-section per CH bond and convolve this with 
a Gaussian profile. These contributions we add. The re- 
sulting envelopes per mode are normalised to the number 
of molecules exhibiting this mode and corrected for the 
number of CH bonds in one functional group. The mea- 
sured cross-sections for the duo modes show a systematic 
decrease with increasing PAH size settling to a value of 
~4 km/mol. To be consistent with the cross-section for the 
largest measured PAHs the contribution of the duo modes 
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Fig. 8. Comparison between the mean interstellar spec- 
trum (panel (b)) and synthetic PAH spectra showing the 
distribution over peak positions of the OOP modes in 



the Hudgins database (Hudgins ct al. 2000b). The shaded 
surfaces in panel (a) and (c) represent the contributions 
per mode for neutral and positively charged PAHs respec- 
tively. Each area represents the average absorption cross- 
section per functional group. 



have been scaled down by a factor 1/4. This way each en- 
velope shows the distribution over wavelengths, while each 
area corresponds to the average absorption cross-section 
per group in the database, see Table |^. Thus the synthetic 
spectrum for each mode is given by: 



Sm{^) — ( G'(A, A„j ^ , Apoak,m,i) 



CH,m 



W„ 



mol.m 



(2) 



where Sm is the synthetic spectrum for mode m, A the 
wavelength, A^^i the cross-section of the mode m in 
molecule i and Apoak.m.i the corresponding peak position. 
G(A,X,r) is a Gauss' function with FWHM=10 cm"!, 
surface=X, peak position=y. A^cH,m designates the num- 
ber of CH bonds in one functional group, A^moLm the total 
number of molecules with mode m. Wm is a weighing fac- 
tor which equals 1, 1/4, 1 and 1 for solo, duo, trio and 
quartet modes respectively to account for the decrease in 
strength of the duo modes in larger PAHs. The summation 
is done over all molecules in the database. The resulting 
envelopes are shown in Fig. ^ for both the neutral and 
cationic PAHs. 



One should bear in mind, that the measured species 
are probably smaller than those that dominate the in- 
terstellar population and that, for stability reasons, the 
interstellar PAH family might be skewed to a few of 
these molecules or the edges structures they represent (cf. 
Sect. |]). Comparing the interstellar spectrum (Fig. ^) 
with these averaged laboratory spectra (Fig. |a,c) allows 
us to further refine the discussion of Sect. 4.1: 



— The position of the strongest band at 11.23 fim agrees 
well with the measured position of solo transitions in 
neutral PAHs but does not agree with the position of 
the cationic solos even after including a 0.1 /im shift 
due to the high temperature of interstellar PAHs. It 
is however important to note that the solo mode of 
the largest cation in the database peaks at 11.2 ^m, 
the longest wavelength of all measured cationic solo 
modes. 

— The peak of the 12.7 /im emission feature falls at 
slightly shorter wavelength than the centre of weight 
for both the neutral and cationic trio modes. Also, the 
blue wing of the 12.7 /im band in the interstellar spec- 
tra does not coincide with any strong emission bands 
of the PAH species in the database. 

— The centre of weight of the neutral quartet vibrations 
is 13.4 /ini which matches better with the position of 
the 13.5 /im UIR band in the interstellar spectrum 
than does the centre of weight for cationic quartets 
(13.25 /im). 

Summarising the above we find that the overall match be- 
tween the UIR spectra and the neutral species is best. 
There are significant differences between the combined 
laboratory measurements and the interstellar spectra. In 
particular the precise assignment of the 12.7 /im feature 
is uncertain. Given the above observations, we feel that 
while the assignment with trio modes is attractive, the 
case in not completely compelling. These issues might be 
resolved when larger species are measured in the labora- 
tory. However it is clear that the 11.0 and 11.2 /im feature 
are due to solo CH bonds in ionised and neutral PAHs 
respectively, and the bands at longer wavelength are due 
to multiplets. 

It is also immediately clear from Fig. || that the in- 
terstellar spectrum does not reflect an equal distribution 
over the different functional groups but is dominated by 
the contribution of solo modes. This reflects the molecular 
structure of the emitting PAHs. 

5. The molecular structures of interstellar PAHs 

It is now possible to quantify the relative amounts of the 
various types of CH groups on the periphery of the inter- 
stellar PAHs which dominate the emission in this wave- 
length range and derive the types and sizes of interstel- 
lar PAH structures implied. This is achieved by analysing 
the peak positions and integrated band strengths listed in 
Table || with the laboratory data summarised in Table |[ 
In the following we assume that the 11.2 /im emission band 
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Table 4. Relative number of solo, duo, trio and quartet 
groups in NGC 7027 and IRAS 18317. 





s/d 


s/t 


s/q 


NGC 7207 


7.7 


4.6 


28.5 


IRAS 18317 


3.4 


0.8 


10.2 



The ratio of the number of solo to duo (s/d), solo to trio (s/t) 
and solo to quartet(s/(7) groups for NGC 7027 and IRAS 18317 
as deduced from their 10 to 15 /im spectra. 




Structure 3. C,2,H35; s/d=6.0; s/t=2.4; s/q=12.0 




Structure 4. CngHj,; s/d=4.0; s/t=1 .1 ; s/q= 8.0 

Fig. 9. Examples of molecular structures simultaneously 

satisfying the structural constraints set by the observed 
band strength ratios of the number of solo, duo and trio 
modes for different interstellar regions. The number of 
solo, duo, trio and quartet functional groups are noted 
s, d, t and q respectively. Solo modes are associated with 
long straight molecular edges. Duos and trios, on the other 
hand, correspond to corners. Quartets are due to pendant 
rings attached to the structure. The numbers in the molec- 
ular structures indicate the number of adjacent CH groups 
per aromatic ring. 



arises from the the solo vibrations, the 12.0 /im band from 
duo modes, the 12.7 /zm from trios and the 13.5 /im from 
quartets. From the laboratory data we get the intrinsic 
strength of solos relative to the duos, trios and quartets 
as 6.2, 2.6 and 2.1 respectively (see Table ||). From the 
intensities of the individual features listed in Table || for 
NGC 7027 we see that the observed solo to duo, solo to 
trio and solo to quartet intensity ratios are 23.8, 4.0, 14.7 
respectively. We conclude that in NGC 7027 the ratio of 
the number of solo CH bonds to CH bonds in duo, trio 



and quartet groups are 3.8, 1.5 and 7.1. Taking into ac- 
count the number of CH's per functional group there are 
7.7, 4.6 and 28.5 solo (groups) relative to the duo, trio and 
quartet groups. For IRAS 18317 these values are (see also 
Table ^) 3.4, 0.8 and 10.2 respectively. These two sources 
represent the extremes in the observed intensity ratio in 
our sample. 

Examples of the types of PAHs which simultaneously 
satisfy these different structural constraints are shown in 
Fig. ^. In constructing these structures one should keep 
in mind that solo H's represent long straight edges, while 
the corners in the structures give rise to duos or trios. 
To match the dominance of the 11.2 /im feature in NGC 
7027 one is naturally driven towards rather large molecules 
with at least ~100— 200 carbon atoms and long straight 
edges. This is entirely in keeping with previous theoreti- 
cal calculations of the molecular sizes of the PAH species 
which account for most of the emission in these features 
( Schutte et al. 1993 ). As illustrated in Fig. ^structure 1 
approximates the ratios listed in Table || for NGC 7027, 
these ratios requires a preponderance of solo hydrogens 
over duos and trios. The extreme quartet to solo group ra- 
tio observed for NGC 7027 is not well reproduced even by 
figure ^ structure 1. In order to reproduce that ratio in one 
single molecule one requires to go to even larger molecules. 
Rather, we surmise the extreme ratio reflects the presence 
of molecules without any quartet groups. Quartet groups 
represent pendant rings on the molecule, which can be 
taken off without altering the other ratios strongly. Of 
course many other PAH structures that match the ob- 
served ratios are possible. However these structures are 
all very similar to this and the conclusion is the same: in 
NGC 7027 the PAH family is dominated by large compact 
PAHs. 

In contrast, for IRAS 18317 the situation is very differ- 
ent. The observed ratios force one to include more corners 
or uneven edges. Structure 4 in Fig. ^ illustrates one way 
of achieving this. Of course this effect can also be achieved 
by going to smaller compact structures of the type shown 
in structure 1 or by breaking up structure 4 in two or more 
fragments. Structures 2 and 3 shown in Fig. ^ are inter- 
mediate between these two extremes and have solo/duo, 
solo/trio, and solo/quartet ratios consistent with the rel- 
ative interstellar band intensities shown in Fig. |^ for the 
objects which lie between the extremes, NGC 7027 and the 
IRAS 18317. Thus we observe a structural evolution where 
closed, compact species dominate the emission in some re- 
gions, while open, uneven structures are more important 
in others. We surmise that this structural evolution as re- 
vealed by the smooth spectral evolution shown in Fig. |^ 
reflects the variations in chemical history and excitation 
environment in these regions. 



6. Discussion 

In this section we will discuss the assignment of 11.2 
and 12.7 interstellar bands and observed variations in 
11.2/12.7 /im ratio and their correlation with other 



12 



Hony et al.: CH OOP bending modes 



UIR bands in view of their underlying physical causes 
(Sect. 



6.1 



In the past these variations have been at- 
tributed to dehydrogenation. We reconsider this sugges 
tion in Sect. p^. 



6.1. The 11.2,12.7 /xm emission features. 

Since quantum calculations show that the CH stretch 
(3.3 /im band) is inherently very weak in ionized PAHs 
(Langhoff 1996), the 3.3 /im emission band is likely a mea- 
sure of the neutral PAH contribution. Based on both the 



peak position of the 11.2 /im band (Sect. 4.2) and its corre- 
lation with the 3.3 /im band (Sect. 3.4) we assign this band 
to solo CH out— of— plane bending vibrations of neutral 
PAHs. This result means that in evolved stars the OOP 
spectrum is dominated by contributions from the neutral 
PAHs. We already assig ned the weak bands at 10.6 and 
11.0 /im to cations (Sect. 4^). Even though the absorption 
cross-sections for the OOP modes do not change upon ion- 
isation, we cannot directly derive the ionisation fraction 
from the measured band strength ratios. This is due to 
the fact that ionisation does strongly affect the strength 
of other modes, which in turn influences the fraction of 



absor bed energy emitted at any wavelength ( Bakes et al, 
2000| ). 



The spectral identification of the 12.7 band is much less 
clear and it is not possible at this time to assign this band 
unambiguously to either neutral or cationic PAHs. In the 
existing database there is only one species with a strong 
band that matches well in position, which 'happens' to be 
a cation. Furthermore the strength of the interstellar 12.7 
/im band correlates with the strength of the 6.2 /im feature 
(see Sect. ^). The strength of the modes between 6 and 
9 /im are greatly enhanced upon ionisation, and thus, one 
way to understand this correlation is to assume that the 
12.7 is also predominantly carried by cations. 

Thus, seemingly, the 11.2 /xm and 12.7 fj,m bands repre- 
sent a dichotomy of interstellar PAHs with the former car- 
ried mainly by neutral and the latter by positively charged 
PAHs. The origin of this interrelation between charge and 
spectral characteristics in unclear. There is no indication 
in the laboratory experiments for a causal relation be- 
tween, for example, charge state and the relative strength 
of the solo to trio modes. Considering also the discussion 
on the molecular structures implied by the relative frac- 
tion of the solos to duos and trios (Sect. Fig. ^ we are 
forced to conclude that the good correlation between the 
11.2 and 3.3 /im bands and between the 12.7 and the 6.2 
/im bands reflects a correlation of molecular structure and 
charge state with environment. Indeed when using PAHs 



containing some 50 C-atoms (Leger & Puget 1984) the 
correlation between the I11.2/I&.2 and 13.3 //e. 2 is well re- 



produced by model calculations of p3akes et al. (2000|) by 
only varying the degree of ionisation. Thus, those envi- 
ronments which favour large PAHs and the 11.2 /im band 
(structure 1 in Fig. ^) also favour neutral PAHs. While in 
regions where open uneven molecular structures and the 



12.7 /im band (structure 4 in Fig. ||) dominate, PAHs are 
predominantly charged. 

This is probably also the origin of the correlation be- 
tween the /pahZ-^ir ratio and the /i2.7/-fii.2 (cf. Fig. ^. 
The IpAn/IiK measures the PAH/dust abundance ratio. 
The loose correlation suggests that for the ISM sources 
the PAH abundance is lower. We recognise that PNe in- 
ject freshly synthesised PAHs into the ISM where they are 
mixed and processed by FUV photons and shocks. This 
processing will lead to a slow destruction of the PAHs. 

The dominant molecular structure reflects the inte- 
grated history of the PAH family and we note that all 
sources with a strong 11.2 /im band are PNe, which 
have formed their PAHs within the last some 1000 years. 
Because open uneven molecular structures are kineti- 
cally more reactive to the addition of carbon atoms than 
compact structures, the predominance of the latter in 
chemically reactive regions where PAHs have recently 



formed can be rationalised. (Frenklach & Feigelson 1989 
phcrchneff et al. 1992 ). In contrast regions with relatively 



strong 12.7 /im bands are all H 11 regions where luminous 
stars illuminate material which has been processed in the 
ISM for some 10^ years. This processing irreversibly leads 
to a breaking down of the molecular structure because 
reformation is prohibited by the low temperature of the 
ISM. 

This does not directly explain why the 11.2 /im 
band correlates with the neutral PAH indicator while 
the 12.7 /im emission feature correlates with bands at- 
tributed to ions. The charge state is rapidly set by the 
charge balance, which is dominated by local physical 
conditions, or more specifically the ionisation parame- 
ter, (Go ^JTgas)/no, where Go is the FUV radiation field, 
Tgas is the gas-temperature and ric the electron density. 
Thus rather than history, ionisation reflects the present. 
Possibly most of the destruction is occurring presently and 
is also driven by local physical conditions. 



6.2. Dehydrogenation 

We have argued that the smooth changes in the band 
strength ratios in the region of interest are caused by 
variations in the edge-structure of the dominant emit- 
ting species. However other effects can also be of influence 
on the emitted spectra. Most notably dehydrogenation. 
There has been a long debate in the literature on the ef- 
fect of dehydrogenation on the spectral characteristics of 
the 10—15 /xm region. Originally when only the 11.2 UIR 
band was known, its dominance had been attributed to ex- 
treme (~90%) dehydrogenation of the emitting aromatic 
species leaving only solo hydrogens (Duley & Williams 



1981). However, this question was revisited when IRAS 



revealed the pres ence of duos and trios in the interstel- 
lar PAH family ( Cohen ct al. 1985 ). Theoretical studies 
of the dehydrogenation of interstellar PAHs have shown 
that for PAHs larger than about 25 C-atoms hydrogena- 
tion through reactions with abundant atomic H is more 
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important than H loss through unimolccular dissociation 
(Tielens et al. 1987; AUamandola ct al. 198£ ; Jochims ct al 



1994 |; |4llain et al. 1996| ; [Jochims et al. 19991 ). Hence, with 
a typical PAH size of 50 C-atoms dehydrogenation should 
have no effect on the UIR spectrum. 

Observationally, our analysis also argues against de- 
hydrogenation. First, we observe a constant ratio of the 
3.3 /im band (all CH oscillation) to the 11.2 /xm band 
(only solo CH oscillation). However, we would expect a 
non-linear behaviour since, when dehydrogenation com- 
mences the number of solo H increases as duos and trios 
are converted to solo's and only at high dehydrogena- 
tion does the relation between the 3.3 and the 11.2 /im 
bands become linear ( ^chutte et al. 1993 ). Secondly, if the 
variation in /12.7//11.2 reflects dehydrogenation than we 
would expect that decreasing H coverage (ie. decreasing 
-'^i2.7/-fii.2) would correlate with increasing CC/CH mode 
emission (ie. /e. 2/^11. 2)- The opposite is actually observed 
(cf. Fig. H). We conclude therefore that dehydrogenation 
has little influence on the observed interstellar UIR spec- 
trum. 



We propose a scenario in which large compact PAHs 
are formed in the winds around evolved stars. These PAHs 
are consequently degraded in the ISM. From the correla- 
tions between charge indicators, which are set by the local 
physical conditions, and the 11.2/12.7 /im band strength 
ratio, which is determined by the molecular structure, we 
conclude that much of this degradation happens on a short 
timescale in the emission objects itself. 
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7. Summary 

We have presented new 10—15 /im spectra of evolved stars, 
H II regions, RNe and YSOs. We observe very rich UIR 
spectra with strong bands at 11.2 and 12.7 /im and weaker 
bands at 10.6, 11.0, 12.0, 13.5 and 14.2 /xm. These spectra 
show large variations in the band strength ratios between 
sources, especially in the 11.2 and 12.7 /xm feature ratio. 
Evolved stars have a dominant 11.2 /im feature while in 
H II regions the 12.7 and 11.2 are typically equally strong. 
We find that the 11.2 /zm band correlates with the CH 
stretch band at 3.3 /im and that the 12.7 /im band corre- 
lates with the CC stretch band at 6.2 /im. 

We have summarised new laboratory spectroscopy re- 
sults for the CH out— of— plane bending vibrations on iso- 
lated neutral and cationic PAHs. Different number of adja- 
cent CH bonds give rise to vibrations in distinctly different 
wavelength regions. The modes are therefore good diag- 
nostics of the molecular structure of the emitting species. 
Upon ionisation the solo CH vibrations are shifted to 
shorter wavelength compared to the solo modes in neu- 
tral. The cross-sections per mode are not strongly modi- 
fied upon ionisation. We attribute the weak bands at 10.6 
and 11.0 /im to solo modes in positively charged PAHs, the 
strong 11.2 /im band the solo modes in neutral. The weak 
12.0 /im band we assign to the duo modes, the 12.7 /im to 
trio modes and the 13.5 /im feature to quartet vibrations. 

From the average cross-sections per mode we have con- 
strained the relative numbers of solo, duo, trio and quartet 
CH groups in different sources for the PAH species that 
effectively emit in this wavelength region. The spectra of 
PNe with a dominant 11.2 /im feature arises from large 
100—150 C-atom) compact PAHs with long straight edges. 
In contrast the H 11 region spectra are due to smaller or 
more irregular PAHs. 
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